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ON ASYMPTOTIC CONSTANTS RELATED TO PRODUCTS OF 
BERNOULLI NUMBERS AND FACTORIALS 



BERND C. KELLNER 



ON 

O 

O, 

^^ . Abstract. We discuss the asymptotic expansions of certain products of Bernoulli num- 

"t^ ' bers and factorials, e.g., 

^; Y[\B2A and n^'"'')''''^ ^^ n^oo 

for integers A: > 1 and r > 0. Our main interest is to determine exact expressions, in 
terms of known constants, for the asymptotic constants of these expansions and to show 
^H ' some relations among them. 

•z.' 

■3 

•^ ' 1. Introduction 

Let Bn be the nth Bernoulli number. These numbers are defined by 
(N 



> • 7 = 2^ Bn-r, Z < 2n, 

e^ — 1 -^^ nl 

n=0 



5£2 ! where -B„ = for odd n > 1. The Riemann zeta function ({s) is defined by 



OO 



)0\ as) = J2^~' = ll^^-P~T\ sec, Res >1. (1.1) 

u=l p 

By Euler's formula we have for even positive integers n that 



an)^-'J^B.. (1,2) 

2 nl 

Products of Bernoulli numbers occur in certain contexts in number theory. For example, 



X 

c^ I the Minkowski-Siegel mass formula states for positive integers n with 8 | n that 



2k J-l 4j, 



which describes the mass of the genus of even unimodular positive definite n x n matrices, 
for details see [12, p. 252]. We introduce the following constants which we shall need further 
on. 
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Lemma 1.1. There exist the constants 

oo 

Ci = JJ((iy) = 2.2948565916... 



u=2 

oo 



C2 = JJ C(2z^) = 1.8210174514... , 

v=l 

00 

C3 = JJ C,{2v + 1) = 1.2602057107... . 

Proof. We have log(l + x) < x for real x > 0. Then 

00 00 00 „ 

lognC(2^) = Y.^oga'^v) < 5^(C(2z.) - 1) = 4- (1.3) 

The last sum of (1.3) is well known and follows by rearranging in geometric series, since 
we have absolute convergence. We then obtain that 7r^/6 < C2 < e^^'^, ({3) < C3 < C2, and 
Ci = C2C3. D 

To compute the infinite products above within a given precision, one can use the following 
arguments. A standard estimate for the partial sum of ({s) is given by 

C(s) - V z/~' < , seR, s>l. 

^-^ s — 1 

This follows by comparing the sum of z/~* and the integral of x~'^ in the interval (A^, 00). 
Now, one can estimate the number A^ depending on s and the needed precision. However, 
we use a computer algebra system, that computes ^(5) to a given precision with already 
accelerated built-in algorithms. Since (^(s) — > 1 monotonically as s — > 00, we next have to 
determine a finite product that suffices the precision. From above, we obtain 

C(s) - 1< 2-M 1 + ^^TT ) , s e M, s > 1. (1.4) 

According to (1.3) and (1.4), we then get an estimate for the remainder of the infinite 
product by 

log n c(^) < 2-^'+^ 

v>N' 

where we can take e = 3/N'; the choice of e follows by 2^ > 1 + xlog2 and (1.4). 

We give the following example where the constant Ci plays an important role; see Finch 
[8]. Let a{n) be the number of non-isomorphic abelian groups of order n. The constant Ci 
equals the average of the numbers a{n) by taking the limit. Thus, we have 

1 ^ 
Ci = lim -— > a(n). 

Af^oo N "^-^ 
n=l 
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By definition the constant C2 is connected with values of the Riemann zeta function on 
the positive real axis. Moreover, this constant is also connected with values of the Dedekind 

eta function 

00 

ry(r) = e'^^'/i^ JJ(1 - e^""'"^), r G C, Imr > 0, 

on the upper imaginary axis. 

Lemma 1.2. The constant C2 is given by 

logp 



I/C2 = JJpi2 rjii- 



P 

where the product runs over all primes. 

Proof. By Lemma 1.1 and the Euler product (1.1) of C{s), we obtain 

00 00 

C2 = n 11(1 -^")"=n 11(1-^")" 

where we can change the order of the products because of absolute convergence. Rewriting 
p~'^'^ = e'^'"^"'^ with r = ilogp /tt yields the result. D 

We used Mathematica [17] to compute all numerical values in this paper. The values 
were checked again by increasing the needed precision to 10 more digits. 

2. Preliminaries 

We use the notation / ~ (? for real-valued functions when lirac^^oo f{x)/g{x) = 1. As 
usual, O(-) denotes Landau's symbol. We write log/ for log(/(x)). 

Definition 2.1. Define the linear function spaces 

fin = spanj^*^, x^'logx}, n > 0, 

0<u<n 

over M where / G f2„ is a function / : IR+ -^ M. Let 

n>0 

Define the linear map ip : f2oo — ^ M which gives the constant term of any / G Qoo- For the 
class of functions 

F{x) = f{x) + 0{x~'), /Gfi„, n>0, 5>0, (2.1) 

define the linear operator [■] : C(]R+;]R) -^ f2oo such that [F] = f and [F] G Qn- Then 
■?/'([-F]) is defined to be the asymptotic constant of F. 

We shall examine functions /i : N ^ M which grow exponentially; in particular these 
functions are represented by certain products. Our problem is to find an asymptotic 
function h : IR+ -^ M where h r^ h. li F = logh satisfies (2.1), then we have [logh] G fin 
for a suitable n and we identify [log h] = [log h] G fin in that case. 
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Lemma 2.2. Let f G Qn where 

n 

f{x) = ^{du x" + Pu x" log x) 

with coefficients a^^Pi, G M. Let g{x) = /(Ax) with a fixed X G M.^ . Then g E Vtn and 
V^(^) = ^(/)+/3ologA. 

Proof. Since g{x) = /(Ax) we obtain 

n 

g{x) = ^(a, (Ax)^ + f3, (Ax)'^(log A + logx)). 

i/=0 

This shows that g G Qn- The constant terms are oq and /^ologA, thus '^/'(f?) = '?A(/) + 
/5ologA. D 

Definition 2.3. For a function / : M.'^ ^ M we introduce the notation 

with functions ]y : M^ ^ M in case / has a divergent series expansion such that 

m— 1 

/(^) = Y. ^-(^) + ^-(^)/-(^)' ^-(^) e (0, 1), m > iV^, 

where iVj is a suitable constant depending on /. 

Next we need some well known facts which we state without proof, cf. [10]. 
Proposition 2.4. Let 

n ^ 

Ho = 0, Hr, = y2-, n>l, 

v=l 

he the nth harmonic number. These numbers satisfy Hn = 7 + logn + 0{n~^) for n > 1, 
where 7 = 0.5772156649... is Euler's constant. 

Proposition 2.5 (Stirling's series). The Gamma function T{x) has the divergent series 
expansion 



logr(x+l) = -log(27r) + (x + - Jlogx-x + J^' 



B2V -(2!^-l) ^ n 

X ^ % X > 0. 



2u{2v- I] 



Remark 2.6. When evaluating the divergent series given above, we have to choose a suitable 
index m such that 



R m-l „ 



^ 2z/(2z/ - 1) ^ 2i.(2i. - 1) + »m{x)n^{x) 

and the remainder \6m{x)Rm{x)\ is as small as possible. Since 6m{x) G (0, 1) is not effec- 
tively computable in general, we have to use \Rm{x)\ instead as an error bound. Schafke 
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and Finsterer [15], among others, showed that the so-called Lindelof error bound L = 1 for 
the estimate L > 9m{x) is best possible for positive real x. 

Proposition 2.7. If a eR with < a < 1, then 

n"" T{n + l-a) 727? fny i_„ 
[U — a) = ;:77^ ^ ~ 777^ T ( — ) 77, 2 ttS U —>■ OO. 



r(l -a) r(l-a) Ve 



Proposition 2.8 (Euler). LetT{x) he the Gamma function. Then 






Proposition 2.9 (Glaisher [9], Kinkelin [11]). Asymptotically, we have 

" 2 

JJi/^ ^ ^ n^"("+^)+i^ e-"^ as n ^ oo 

where A = 1.2824271291... is the Glaisher-Kinkelin constant, which is given by 
log.4 = l-C'(-l) = ^4log(2.)-^. 

Numerous digits of the decimal expansion of the Glaisher-Kinkelin constant A are 
recorded as sequence A074962 in OEIS [16]. 

3. Products of factorials 

In this section we consider products of factorials and determine their asymptotic expan- 
sions and constants. For these asymptotic constants we derive a divergent series represen- 
tation as well as a closed formula. 



Theorem 3.1. Let k be a positive integer. Asymptotically, we have 

|n(n+l) 

CX) 

u=l 

with certain constants Tk which satisfy 



W{ku)\ ^ TkA\2'K)-*i^\ (27rA:e'=/2-i^^t^i+A+_i_ ^^^ 



WF - ^ I V- B^^ ^(2j - 1) 
^ "^ 12fc 4^ 2j(2j-l)F^-i' 

Moreover, the constants have the asymptotic behavior that 

n 

lim jFjt = 1, lim JF^ = e"''^'^ ^ and 1 I -^fe ~ ^oo rC '^ as n ^ oo 

with 



fc— >oo fc— »oo 

k=l 



log-^c« = y^ + 5^ 



Y , V-, B^, C(2j - 1 



12 ^ 2,(2, 
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Theorem 3.2. If k is a positive integer, then 

We will prove Theorem 3.2 later, since we shall need several preparations. 

Proof of Theorem 3.1. Let fc > 1 be fixed. By Stirling's approximation, see Proposition 
2.5, we have 

\og{ku)\ = ^ log(27r) + (ku + 1") \og{kv) - ku + f{ku) (3.1) 



where we can write the remaining divergent sum as 



B2j 



12ku ^ 2j(2j-l)(H2^-i- 



Define S{n) = 1 + ■ ■ ■ + n = n{n + l)/2. By summation we obtain 

y \og{ki')\ = — log(27r/c) H — logn! — kS{n) + kS{n) log k 

n n 

k^ulogu + ^f{ku). 



2 



The term |logn! is evaluated again by (3.1). Proposition 2.9 provides that 

k 
V log u = k log A + kS[n) log ri + 



i/logz/ = A;log^+ kS{n) logn -\ logn (5'(n) j + 0{n~ 



with some 5 > 0. Since lim^^oo -f^n — log ra = 7, we asymptotically obtain for the remaining 
sum that 



n 



^2k ^ 12k ^ 



hm y f{ku) -^\ogn]=^ + y' T '^ ^T;/ , =: log^^. (3.2) 

r^oc\^''^ ' I2k ^ / 12A; ^^ 2? 27 - lU^J-i ^ ^ ^ 



2j(2j 



Here we have used the following arguments. We choose a fixed index m > 2 for the 
remainder of the divergent sum. Then 



r V^fl (h \ -"2m _ -P2mU2m- ij 

i™ Z.'^-^'^'^^2m(2m-l)(H''"-' " ^-2m(2m - 1) P— 1 ^ ^ 
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with some rj^ G (0, 1), since 9m{kv) G (0, 1) for all z/ > 1. Thus, we can write (3.2) as an 
asymptotic series again. Collecting all terms, we finally get the asymptotic formula 

" 1 / 3 \ 

^\og{kv)\ =logJ^fc + A:log^+-log(27r)+ kS{n) ( -- + log(A;n) J 

+ ^nog(27rA;) + --l+logn 

- H \ r logn + Oin ^ ) 

with some 6' > 0. Note that the exact value of 6' does not play a role here. Now, let k be 
an arbitrary positive integer. From (3.2) we deduce that 

logTk = ^ + 0{k-^) and klogJ^k = Y^ + 0{k~'). (3.4) 

The summation of (3.2) yields 



Eiog^^^i/,.EE' ,:;^^r;;.i. (3.5) 

fc=l fc=l j>2 



Similar to (3.2) and (3.3), we can write again: 

12--'^ ) 12 ^ 2j(2j 



■rjE'-^-^'o.") -^-E'^|i^-o^^». (3.«^ 



n— »oo 



The case /c = 1 of Theorem 3.1 is related to the so-called Barnes G-function, cf. [2]. 
Now we shall determine exact expressions for the constants JF^. For k > 2 this is more 
complicated. 

Lemma 3.3. We have Tx = (27r)i eT2/A'^. 

Proof. Writing down the product of n\ repeatedly in n + 1 rows, one observes by counting 
in rows and columns that 

n n 

n!"+i = JJ z/! JJ z/^ (3.7) 

From Proposition 2.5 we have 

n + 1 f 1\ 1 
{n + l)logn! = log(27r) -n{n + 1) + (n + 1) ( n + - j lognH V0{n~^). 

Comparing the asymptotic constants of both sides of (3.7) when n ^ cxd, we obtain 

{2TT)2e^ = J^^A{2Ti)-i -A 
where the right side follows by Theorem 3.1 and Proposition 2.9. D 
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Proposition 3.4. Let k,l be integers with k>l. Define 



Fk,i{n):=f[{kiy-l)\ for < I < k. 



u=l 



Then [logF^,/] G ^2 and Fkfl{n) ■ ■ ■ Fk,k-i{n) = Fifi{kn). Moreover 



FkAn)lFk,i+i{n) = k^\[ 



v=l 



z/--) for Q<1 <k-l 



and [\og{Fk,i/Fk,i+i)] = [log Fk,i] - [logFfc,/+i] G l^i. 

Proof We deduce the proposed products from {ku — l)\/{ku — (/ + 1))! = ku — I and 



kn 



W{kv)\{kv - 1)! ■ ■ ■ {kv -{k- 1))! = JJ v\. 



(3J 



u=l 



u=l 



Proposition 2.7 shows that [Iog(Ffc //F^ /+i)] G fli. Since the operator [■] is linear, it follows 
that 

[\og{Fk,i/Fk,i+i)] = [log Fk,i - log Fk,i+i] = [log Fk,i] - [log Fk,i+i] G Qi. (3.9) 

From Theorem 3.1 we have [logF^o] G ^2- By induction on / and using (3.9) we derive 
that [logFfc,;] G (^2 for < / < A;. ' D 

Lemma 3.5. Let k be an integer with k > 2. Define the k x k matrix 

\ 



M, : = 



/I -1 

1 -1 



\1 1 



1 -1 
1 1 / 



where all other entries are zero. Then det Mk = k and the matrix inverse is given by 



M-^ = \Mk with 



( k-\ k-2 k-3 
-1 k-2 k-3 
-1 -2 k-3 



Mk = 



2 
2 
2 



1\ 
1 

1 



V 



-2 -3 
-2 -3 
-2 -3 



2 1 1 

-{k-2) 1 1 

-(fc-2) -{k-1) 1 J 



Proof. We have det M2 = 2. Let A; > 3. We recursively deduce by the Laplacian determi- 
nant expansion by minors on the first column that 

det Mk = (-1)^+^ det Mfc„i + (-1)^+'= det Tk-i 
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where the latter matrix T^-i is a lower triangular matrix having —1 in its diagonal. There- 
fore 

det Mfc = det Mfc_i + (-1)^+^' ■ (-1)'="^ = k-l + l = k 

by induction on k. Let /^ be the k x k identity matrix. The equation M^ ■ M^ = k I^ is 
easily verified by direct calculation, since Mk has a simple form. D 

Proof of Theorem 3.2. The case k = 1 agrees with Lemma 3.3. For now, let k > 2. We 
use the relations between the functions F^/, resp. logF^;, given in Proposition 3.4. Since 
[logFfc/] e ^2; we can work in ^2- The matrix M^ defined in Lemma 3.5 mainly describes 
the relations given in (3.8) and (3.9). Furthermore we can reduce our equations to M by 
applying the linear map ip, since we are only interested in the asymptotic constants. We 
obtain the linear system of equations 



where 



M.-x = b, x,6gM*^, 



x = (^([logFfc,o]),...,^([logFfc,,_i]))^ 
and b = (6i, . . . , bk)^ with 

6,+i = ^P{[\og{Fk,i/Fk,i+i)]) = ^ log(27r) - log F Tl - ^"j for l = 0,...,k-2 

using Proposition 2.7. The last element bk is given by Theorem 3.1, Lemma 3.3, and 
Lemma 2.2: 

bk = ^([log(Fi,o(A;n))]) = - log(27r) + log^i + log^ + — logk 
= - log(27r) - log ^ + — + — log k. 



By Lemma 3.5 we can solve the linear system directly with 
The first row yields 



x = -Mk-b. 



''^ = T.^>^ + T.H^^-'''>^'^ 



v=l 



On the other side, we have 



This provides 



xi = ^([logFfc,o]) = logJ^fe + - log(27r) + fclog A 



log J-, = - ( A: + 1) log^ + (^^ + ^ - i ) log(27r) 

fe-i 



1 v;^ V — \ (V 



Vlk 

u=2 



(3.10) 
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after some rearranging of terms. By Euler's formula, see Proposition 2.8, we have 

^i:i°^r(9 = (i-±)log(2.)-llog.. (3.11) 

u=l ^ ' 

Finally, substituting (3.11) into (3.10) yields the result. D 

Remark 3.6. Although the formula for JF^ has an elegant short form, one might also use 
(3.10) instead, since this formula omits the value r{l/k). Thus we easily obtain the value 
of JF2 from (3.10) at once: JF2 = (27r)4 224 e^s /A^ . 



Corollary 3.7. Asymptotically, we have 

with the constants e^/A = 0.8077340270... and {2tt)^/A = 1.2345601953.... 
Proof. On the one hand, we have by (3.4) that 

nlogJ^n = T^ + 0{n~^). 
On the other hand. Theorem 3.2 provides that 

T 1 2 "-1 

nXogTn = - (n^ + l)log^+ — - — logn + — log(27r) - Vi/logr(- 

Combining both formulas easily gives the result. D 

Since we have derived exact expressions for the constants JF^, we can improve the calcu- 
lation of Too- The divergent sum of JF^o, given in Theorem 3.1, is not suitable to determine 
a value within a given precision, but we can use this sum in a modified way. Note that we 
cannot use the limit formula 

/ n. 

7 



logJF^ = lim V'logJF^ - — logn 

77, — >r)C) \ ' * /. 



n^oc \ ' ' 12 

.fc=l 



without a very extensive calculation, because the sequence 7„ = iJ„, — log n converges too 
slowly. Moreover, the computation of T}^ involves the computation of the values Y{vjk\ 
This becomes more difficult for larger k. 

Proposition 3.8. Let m, n be positive integers. Assume that m > 2 and the constants J^k 
are given by exact expressions for k = 1, . . . ,n. Define the computable values r]k € (0, 1) 
implicitly by 

, 7 ^' i?2jC(2j-l) i^2mC(2m-l) 

°^ *^ 12fc ^ 2j{2j - 1) Fi-i '''^2m(2m - 1) F--^' 
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Then 

2 m-l 



iogy-oo = 77: 



12 ^ 2j{2j-l) "'"" 2m(2m-l) 

w^th ^„,„ G (^^j^,^^^^) C (0, 1) where 

0- = C(2m - 1)- J2 ^, 9- = 1 + C(2m - 1)- J] |^. 

fc=l A:=l 

T/ie error bound for the remainder of the divergent sum of log J-'oo is given by 

152^1 C(2m- 1)2 



On!m = 1 - C(2m 



k=l / 



2m(2m 



Proof. Let ri > 1 and m > 2 be fixed integers. The divergent sums for log^^ and log^oo 
are given by Tlieorem 3.1. Since we require exact expressions for JF^, we can compute tlie 
values r]k ioT k = 1, . . . ,n. We define 

Vm,k = v'm,k = Vk for A; = 1, . . . , n 
and 

'7m,fc = 0, r/'^ ;, = 1 for k > n. 
We use (3.5) and (3.6) to derive tlie bounds: 



OO OO / 



J^2m-l "''^" ^V-"" -/ _^;;_^ ^2m-1 "n,m- 

fe=l fc=l 



We obtain the suggested formulas for ^™™ and ^™^ by evaluating the sums with r]m,k = 0, 
resp. rj'^i^ = 1, ioY k > n. The error bound is given by the difference of the absolute values 
of the minimal and maximal remainder. Therefore 



^ Z_^ j/2m-l j 
k=l / 



withi?= |52^|C(2m-l)V2m(2m-l). D 

Result 3.9. Exact expressions for J^j.: 

15 1 



J^l = (27r)4eT5/^2^ J^2 = {271)^2^6^ /A 

,5 5 1,,10 /n\i , ,1 1 1,,17 /ox 

J^^ = (27^)T5 3^e^/^-^(|)^ J^4 = (27r)2 23 e^/^x r(f) 



We have computed the constants J-'k by their exact expression. Moreover, we have 
determined the index m of the smallest remainder of their asymptotic divergent series and 
the resulting error bound given by Theorem 3.1. 
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Constant 


Value 


m 


Error bound 


^1 


1.04633506677050318098... 


4 


6.000 ■ 10-^ 


T2 


1.02393741163711840157... 


7 


7.826 ■ 10-^ 


^3 


1.01604053706462099128... 


10 


1.198-10-9 


T, 


1.01204589802394464624... 


13 


1.948 ■ 10-12 


^5 


1.00963997283647705086... 


16 


3.272 ■ 10-1^ 


-^6 


1 .00803362724207326544. . . 


20 


5.552 ■ 10-18 



The weak interval of Too is given by Theorem 3.1. The second value is derived by 
Proposition 3.8 with parameters ?7i = 17 and n = 7. Thus, exact expressions of T\^ . . . ,Tj 
are needed to compute Too within the given precision. 



Constant 


Value / Interval 


m 


Error bound 


''00 
•J 00 


(1.02428, 1.02491) 
1.02460688265559721480... 


4 
17 


6.050 ■ 10-^ 
6.321 ■ 10-22 



4. Products of Bernoulli numbers 

Using results of the previous sections, we are now able to consider several products of 
Bernoulli numbers and to derive their asymptotic expansions and constants. 



Theorem 4.1. Asymptotically, we have 

n \ n{n+l) 



n 1^2.1 


~ Bi 


Y' -^21/ 


~ B2 


11 2u 


with the constants 





7j-g3/2 

n 



(167m) 2 77.24 as n 



CXD, 



vre 



3/2 



An\ 2 
Tie J 



n 24 as n 



00 



Bi = C2J^2^2(2vr)3 = C2 (27r)l 2 A e^/^i, 
B2 = C2J'2^2/(27r)i = C2 2li e^Mi 

Proof. By Euler's formula (1.2) for C(2z^) and Lemma 1.1 we obtain 

2-f2z/)! 



n 1^2.1 ~ C2I1 



J- (27r) 
Theorem 3.1 states for k = 2 that 

n 



2v 



C2 2"(27r) 



-n(n+l) 



11(2^)' 



as n 



00. 



v=\ 



J](2z/)! ~ T2A\2tx)\ 



v=\ 



2n 



n{n+l) 



(47rn) 2 n24 as n 



00. 



The expression for T2 is given in Remark 3.6. Combining both asymptotic formulas above 
gives the first suggested formula. It remains to evaluate the following product: 



n(2 

i/=i 



2"n\ ~ (27r)5/'— Vn^ 



as n 



CXD. 
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After some rearranging of terms we then obtain the second suggested formula. D 

Remark 4.2. Milnor and HusemoUer [14, pp. 49-50] give the following asymptotic formula 
without proof: 

n 

Y[\B2u\ ~ i3'n!2"+^F(2r2 + l) as n ^ oo (4.1) 

where 

2 / \ — 

and B' ~ 0.705 is a certain constant. This constant is related to the constant S2. 

Proposition 4.3. The eonstant B' is given by 

B' = 2^ 2-i B2 = C2 e^/2^A^ = 0.7048648734... . 

Proof. By Theorem 4.1 we have 

'"' \R \ 
JJ^M ^ B2G{n) as n ^ oo (4.3) 



2u 



with 



2 / 1 



We observe that (4.1) and (4.3) are equivalent so that 

2B'F{2n + l) ~ B2G{n) as n -^ oo. 
We rewrite (4.2) in the suitable form 



F(2„+i) = (";f' 


"'+'^+ V 47re \ t+i / . / 1\^ 

/ O o,i /n 1 




Hence, we easily deduce that 




G(n)/F(2n+1) = 


IH 62 224 

\ 2n) V 4 7 


It is well known that 




lim i\ +-V = P^ 


/ X\"^ 2!2 

and lim 6"^'" M + - = 6"~ 



n— >oo \ 11/ n^xxi \ fl/ 

Evaluating the asymptotic terms, we get 



1 111 



2B'/B2 ~ G{n)/F{2n + 1) ~ 68 6-4 224 68 2-2 as n ^ 00, 
which finally yields i3' = 2242"2i52. □ 
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Theorem 4.4. The Minkowski-Siegel mass formula asymptotically states for positive in- 
tegers n with 4 I n that 



M(2n) 



2u\ 



I R I "~^ I R 

2n J-l 4i, 

i/=i 



i3. 



n 



ne 



3/2 



4ra\ 2 ^ 

— n24 as 72 — > oo 

Trey 



Proof. Let n always be even. By Proposition 2.5 and (1.2) we have 



Bn/n 



B2n/2n 
since C('^)/C(2^) ~ 1 ^^^ 

, , n! 
log 



O C(^) ^r. "-' 
2 — : -TT 



(2n)! 



C(2n) (2n)! 



n — n log n 



^(S) 



as n ^ oo, 



2n H — ) log 2 as n ^ cxd. 



We finally use Theorem 4.1 and (4.3) to obtain 



M{2n) 



which gives the result. 



Bn/n 



B2n/2n 



n 



u=l 



\B. 



2u\ 



2v 



\7re 



G{n) as n -^ oo, 



D 



Result 4.5. The constants B', By [v = 1, 2, 3) mainly depend on the constant C2 and the 
Glaisher-Kinkelin constant A. 



Constant 


Expression 


Value 


A 




1.28242712910062263687... 


C2 




1.82101745149929239040... 


Bi 


C2(27r)5 2^6 21/^5 


4.85509664652226751252... 


B2 


C22A 6 21/^3 


1.93690332773294192068... 


B, 


^ 17 1 / . 1 
C2224 6 24/^2 


2.73919495508550621998... 


B' 


C2 623/2!^^ 


0.70486487346802031057... 



5. Generalizations 

In this section we derive a generalization of Theorem 3.1. The results show the structure 
of the constants J-'k and the generalized constants J>,fc, which we shall define later, in a 
wider context. For simplification we introduce the following definitions which arise from 
the Euler-Maclaurin summation formula. 

The sum of consecutive integer powers is given by the well known formula 



n-l 

E 

i/=0 



u 



Br+i{n) - Br+l 
r + 1 



SG)^-7^' 



r > 0, 
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where B^^x) is the ??7,th Bernoulh polynoiniaL Now, the Bernoulh number i?i = — | is 
responsible for omitting the last power n^ in the summation above. Because we further 
need the summation up to n^, we change the sign of Bi in the sum as follows: 

^.(n) = 5^.^ = 5^r (-ir-i?._,^, r>0. 



■,„+! -Dn+1 



This modification also coincides with 

n + 1 
for nonnegative integers n. We define the extended sum 

where the symbol o is replaced by the index j ' + 1 in the sum. Note that 5*,. is linear in the 
second parameter, i.e., 

Sr{n; a + /3/(o)) = aSrin) + /5S',(n; /(o)). 

Finally we define 

Dkix) = y^' IBon kX~^^-'~^' where B^k = — 7 —n r- 

^ ^ ^ ^' ' m(m - l)k"'-^ 

Theorem 5.1. Let r be a nonnegative integer. Then 

n 

llz/'' ~ ArQr{n) as n ^ oo. 



u=l 



where Ar is the generalized Glaisher-Kinkelin constant defined by 

\ogAr = -a-r)Hr-C{-r). 

Moreover, logQr G ^r+i with 

\ogQr{n) = {Sr{n) - C{-r)) logn + Sr{n] H,. - H^). 

Proof. This formula and the constants easily follow from a more general formula for real 
r > — 1 given in [ID, 9.28, p. 595] and after some rearranging of terms. D 



Remark 5.2. The case r = reduces to Stirling's approximation of n\ with ^o = v27r. 
The case r = 1 gives the usual Glaisher-Kinkelin constant Ai = A. The expression 
Sr{n; Hr — H^) does not depend on the definition of Bi, since the term with Bi is cancelled 
in the sum. Graham, Knuth, and Patashnik [10, 9.28, p. 595] notice that the constant 
— (^'(— r) has been determined in a book of de Bruijn [7, §3.7] in 1970. The theorem above 
has a long history. In 1894 Alexeiewsky [.3] gave the identity 

n 

n^^^=exp(C'(-^,n + l)-C'(-r)) 
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where C'{s, a) is the partial derivative of the Hurwitz zeta function with respect to the 
first variable. Between 1903 and 1913, Ramanujan recorded in his notebooks [5, Entry 27, 
pp. 273-276] (the first part was published and edited by Berndt [5] in 1985) an asymptotic 
expansion for real r > — 1 and an analytic expression for the constant Cr = —('{—r). 
However, Ramanujan only derived closed expressions for Co and C2r (r > 1) in terms of 
({2r + 1); see (5.9) below. In 1933 Bendersky [4] showed that there exist certain con- 
stants Ar- Since 1980, several others have investigated the asymptotic formula, including 
MacLeod [13], Choudhury [G], and Adamchik [1, 2]. 

Theorem 5.3. Let k, r be integers with k > 1 and r > 0. Then 

n 



u=l 



The constants J^r,k o.iT'd functions Pr^k satisfy that lim J^r,k = 1 o.nd log P^. ^ G ^2^+2 where 

fc— >oo 

\ogPr,k{n) =-Sr{n)\og{27ik) + k Sr+i{n)\og{k/e) 

+ Br+2,k logn + ^ B2j,k Sr+l-2j{n). 

i=i 
The constants Ar and functions Qr are defined as in Theorem 5.1. 

The determination of exact expressions for the constants J^r,k seems to be a very com- 
plicated and extensive task in the case r > 0. The next theorem gives a partial result for 
k = 1 and r > 0. 

Theorem 5.4. Let r he a nonnegative integer. Then 



\0gTr,l = 2 log A - log A'+l + Sr{l; fil+o,l " log A 



Case r = 0: 
Case r > 0: 

where 



logj; 



r,l 



12 



logA-21ogA 



r+l 



a 



r,J 



log J^,. 1 = arfi + ^ arj logAj 
i=i 

2!^, r^j (mod 2), j = 0; 

tOMm^ r.,(mod2), , = 0; 
3=0 



-s.. 



.+!,,- (7)^^, r^j (mod 2), j>0; 
0, r = j (mod 2), j > 



and 6ij is Kronecker's delta. 
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Proof of Theorem 5.3. Let k and r be fixed. We extend the proof of Theorem 3.1. From 
(3.1) we have 

\og{ku)\ = -\og{27rk) + kulog (-] + (ku + -] logu + Dkiu). (5.1) 



The summation yields 



where 



Y^ ^" \og{kv)\ = Fi{n) + F2(n) + F^{n) 



Fi{n) = -Sr{n)\og{2nk) + kSr+i{n)\og{k/e) 



n 1 '^ 

F2{n) = kY^ v^^^ logv + '^v'' log , 

n 

Fs{n) = J2^^M^)- 



u=l 



Theorem 5.1 provides 



F2{n) = k (log A+i + logQr+i{n)) + - (log A + logQrin)) + 0{n~^) 



with some 6 > 0. Let R = [-^J- By definition we have 



2 
R 



j=l j>R 

Therewith we obtain that 

R n 

^3('^) = X] ^2i,fc Sr+l-2jin) + ^ E2{u). 
j=l u=l 

For the second sum above we consider two cases. We use similar arguments which we have 
applied to (3.2) and (3.3). If r is odd, then 

n 

hm J2e^(^) = Z'^2,>C(2j - (r + 1)). (5.2) 

n— >oo ^— ' •^— ' 

u=l j>R 

Note that Br+2,k = in that case. If r is even, then we have to take care of the term v~^. 
This gives 

lim (y^E2{u)-%^2,k\ogn\ =7^r+2,fc+ V'52,,fcC(2j - (r + 1)). (5.3) 

\i'=l / j>R+l 
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The right hand side of (5.2), resp. (5.3), defines the constant logj-'r^k- Finally we have to 
collect all results for Fi, F2, and F^. This gives the constants and the function Pr,fc- It 
remains to show that linifc_^oologJ^r,fe = 0. This follows by B2j,k — > as /c ^ 00. D 

The following lemma gives a generalization of Equation (3.7) in Lemma 3.3. After that 
we can give a proof of Theorem 5.4. 

Lemma 5.5. Let n,r be integers with n > 1 and r > 0. Then 



n\ 

v=l u=l v=l 



2^' 




31^ 


22" 




32' 



Proof. We regard the following enumeration scheme which can be easily extended to n 
rows and n columns: 

12'- 

13'- 2^^ 

The product of all elements, resp. non-framed elements, in the z/th row equals n!'^'^, resp. 
ul^"^ . The product of the framed elements in the z/th column equals u^^^'^~^\ Thus 

n n 

v=l v=l 

Proof of Theorem 5.4- Let r > 0. We take the logarithm of (5.4) to obtain 

Fi(n) + F^in) = F^in) + F^{n) (5.5) 

where 

n 

Fi{n) = Sr{n)\ogn\, F2{n) = \Jz/''logi 



v=l 
n 

v=l v=l 



F3(n) = J^z/Hogz/!, F^{n) = Y,SrWog^ 



Next we consider the asymptotic expansions Fj of the functions Fj (j = 1, ... ,4) when 
n — i> 00. We further reduce the functions Fj via the maps 

C(M+;M) Al^ Q^ JU R 

to the constant terms which are the asymptotic constants of [Fj] in fioo. Consequently 
(5.5) turns into 

^pm]) + ^pi[F2]) = iJim) + ij{[F,]). (5.6) 

We know from Theorem 5.1 and Theorem 5.3 that 

^Pi[F2]) = log A and ^^([Fg]) = log J^,,i + - log A + log A+i- 
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For F4 we derive the expression 

ij{[F,]) = Sr{l;logA,), (5.7) 

since each term SjU^ in ^^(z/) produces the term SjlogAj. It remains to evaluate Fi. 
According to (5.1) we have 

logra! = - log(27r) — n + in -\ — J logn + Di{n) =: E{n) + Di{n). 

Thus 

Fi(x) = Sr{x)E{x) + Sr{x)Di{x). 

Since S^E G fioo has no constant term, we deduce that 

^([Fi])=V^([5,Di]) = ^,(l;fii+,,i). 

The latter equation is similarly derived as (5.7), whereas we regard the constant terms of 
the product of the polynomial Sr and the Laurent series Di. From (5.6) we finally obtain 

log j;,i = - log A - log A'+i + Sr{l] 5i+o,i - log A). 

Now, we shall evaluate the expression above. For r = we get 

1 1 
logJ^o.i = Y^ + 2 ^°^^° - 2 log A, 

since 

5*0(1; 5i+o,i - log A) = ^2,1 - log A = T^~ log^i- 

For now, let r > 0. We may represent logJ?>i in terms of log A ^^ follows: 

r+l 
logjF^^l = ar,o + ^^Olr^j log A- 
j=l 

The term a^fl is given by 

a.,0 = 5.(1; fii+.,i) = Y. {"^ (-i)-^fi,_^.:^ 

j=o ^-^^ ^ 

where the sum runs over even j, since -Bj+2,1 = for odd j. If r is odd, then the sum 
simplifies to the term i?j.+i/2r(r + 1). Otherwise we derive for even r that 



a, 



El r\ Br^j Bj^2 



,=0 j7(j + iP(j + 2)' 



It remains to determine the coefficients arj for r + 1 > j > 1 . Since ^x^ — x^~^^ — Sr (x) is an 
odd, resp. even, polynomial for even, resp. odd, r > 0, this property transfers in a similar 
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way to i log Ar — logAr+i — Sr{l] log^o), such that a^j = when 2 | r — j. Otherwise we 
get 

f r \ Br-(j-i) [r + l\ Br+i^j ^ ,_ ^, 

''r, = -[^_^)^f^-Sr+l,j = -(^ ^ )^:^-Sr+l,j (5.8) 

for 2 \ r — j, where the term — log^^+i is represented by —6r+ij- □ 

Corollary 5.6. Let r be an odd positive integer. Then 



1 — 1 

2 



.0.-....^-^ ^4c(-i-.)c(.^i,-^^:±^f:±^ 



1 — 1 

2 



-1) — ^ 



ir! / \Br+i\ ^ \Br+i-2j\Ci'2j + l) (r + 2)C(r + 2 ) 
2 1 r(r + 1)! ^ (r + 1 - 2j)\ (27r)2i (27r)'-+i 



Proof. As a consequence of the functional equation of ({s) and its derivative, we have for 
even positive integers n, cf. [5, p. 276], that 

1 n' 
logX = -C'(-n) = -^J^-y^Cin + 1) (5.9) 

where the left hand side of (5.9) follows by definition. Theorem 5.4 provides 

r+l 

log j;,l = /^\^ + y]ar,2j log^2j. 

2r r + 1) -^^^ 

Combining (5.8) and (5.9) gives the second equation above. By Euler's formula (1.2) we 
finally derive the first equation. D 

Remark 5.7. For the sake of completeness, we give an analogue of (5.9) for odd integers. 
From the logarithmic derivatives of r(s) and the functional equation of C(s), see [5, pp. 183, 
276], it follows for even positive integers n, that 

log A-i = — i/n-i - C'(l -n) = ^(7 + log(2vr)) + 2^-^^-^C{n) 
where 

oo 



C'('^) = -^^og{iy)iy- 



u=2 



However, Mathematica is able to compute values of C for positive and negative argument 
values to any given precision. 

Result 5.8. Exact expressions for jF^i in terms of ^j: 
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Constant 


Expression 


Value 


-^0,1 

-^1,1 

•^3,1 


7 -1 ** 
6540^^ M3 ' 

1 _1 _5 

g7 1 1 6 
e 18900^f ^g^^gS 

6 2520 ^2'' A ''A' 


1.04633506677050318098... 
0.99600199446870605433... 
0.99904614418135586848... 
1.00097924030236153773... 
1.00007169725554110099... 
0.99937792615674804266... 



Exact expressions for Tr,\ in terms of C(2j + 1): 

' 1 3C(3)\ 
24 87r2 )' 



exp 



J^3,i = exp 

JFsi = exp 
For the first 15 constants Tr 1 (r 



1 
720 
1 



C(3) 
167r2 

C(3) 

2520 ^ 487r2 

= 0, . . . , 14) we 

11 



max \Tr\ - 

0<r<14 

but, e.g., ^19^1 K, 371.61 and J^2o,i ~ 1-16 ■ 10 



15C(5) \ 
167r4 y' 

505) _ 10507) 
167r4 167r6 

find that 

< 0.05, 
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